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Abstract 

We verify whether the previously developed model of a KD2PO4 crystal, with the shear strain 
£6 taken into account, is able to describe the longitudinal electric field influence on the KH2PO4 
family ferroelectrics. Major effects of the strain eq are splitting of the Slater energies of the short- 
range correlations and the effective field created by piezoelectric coupling. Calculated Tc — -E3 phase 
diagrams, field dependences of polarization, susceptibility, and elastic constant of deuterated KD2PO4 
and pure KH2PO4 well accord with the available experimental data. For a consistent description of 
all dielectric and piezoelectric characteristics of the crystals, phonon degrees of freedom must be taken 
into account. 

1 Introduction 

Studies of ferroelectric crystals behavior under different external factors such as pressure or electric field 
provide an important information about the phase transition mechanisms in these crystals as well as of 
their dielectric and piezoelectric responses. Thus, studies of the transition temperature dependence on 
the geometry of a hydrogen bond revealed a universal dependence of the transition temperature on the 
distance between equilibrium proton (deuteron) positions on a bond in several crystals of the KH2PO4 
family ^j. Studies of the effects produced by external fields conjugate to the order parameter (shear 
stress (Tg a-nd electric field E3 in these crystals) permit to elucidate the role of structural changes occurring 
at the phase transition, such as appearance of the shear strain £g. 

In our recent paper |^ we modified the proton ordering model for deuterated KD2PO4 type crystals 
so that the shear strain eg was taken into account. We obtained a fair agreement of the theoretical results 
for associated with strain eg dielectric, elastic, and piezoelectric characteristics of KD2PO4 (piezoelectric 
constants d^e, 635, /i3g, (/ag, elastic constants Cqq, cfg and compliances Sqq, sfg, and dielectric permittivity 
£33) with the relevant experimental data. 

A consistent taking into account of the changes caused by the strain eg in the system symmetry 
revealed that this strain splits the energies of short-range deuteron configurations. Thus, the lowest level 
of the so-called up and down configurations and the energies of lateral and single-ionized configurations 
are split. Splitting of up/down and ionized configurations gives rise to single- and three-particle spin terms 
in the short-range Hamiltonian. Our numerical calculations have shown that by fitting to experimental 
data one is not able to determine the magnitude of the splitting of up/down configurations ^^g, since if a 
certain ratio between S^q and the magnitude of piezoelectric internal field -^g Q , the experiment is well 
described in a rather wide range of dss- A good agreement with the experimental data for temperature 
behavior of the above mentioned physical characteristics of KD2PO4 is also obtained in a simplified 
version of the model ||^, where we neglected splitting of the up/down and ionized configurations and took 
into account only splitting of lateral configurations and piezoelectric internal fields. 

In we show that the transition temperature of KD2PO4 type crystals is increased by the stress erg, 
and the jumps of the order parameter are decreased. At the constructed phase diagram there are two 



critical points where the phase equilibrium curves terminate. Stresses above critical smear off the phase 
transition and lead to smooth temperature dependences of polarization and strain. Correspondingly, the 
peak values of those characteristics of the crystal that have peculiarities at the transition points (the 
longitudinal dielectric permittivity, compliance s^g, piezomodules d^Q and 630, and specific heat) increase 
with stress. These peaks are the highest at the critical stress, whereas at stresses above critical their 
temperature curves are smooth. Such a behavior is peculiar to the ferroelectrics with the first order phase 
transition in external fields conjugate to the order parameter. 

Studies of the stress ae influence on the KH2PO4 family crystals, though interesting they are, remain 
quite abstract, since their experimental verification appears to be very difficult. Instead, measurements 
in electric fields are a more accessible experimental technique. The longitudinal electric field E3 should 
be considered, since its influence on the physics of the phase transition the relevant characteristics of the 
KH2PO4 type crystals is analogous to the influence of the shear stress ae and is well studied experimentally 
[101,1,0,0,11,0. 

In this paper we shall show that the developed theory is capable, without introducing any additional 
fitting parameter, of describing the behavior of the phase transition, polarization, and other characteristics 
of a deuterated KD2PO4 in electric field £'3, and attempt to apply this theory also for the description of 
the relevant characteristics of a pure KII2PO4, with tunneling neglected. 



2 Phenomenological analysis 

Strictly speaking, for the KH2PO4 type crystals a expansion of thermodynamic potential G in polarization 
P3 and strain eg 



G = Go + - {aiPi + a2P3£6 + aseg) + 7 (&i^3 + ^2-P|^6 + 63P3 ^e) 



+ i {ciPl + C2Piee + csP^el) - VE3P3 - va^e^ (2.1) 

is the most consistent. It allows one to perform the analysis for the both fields conjugate to the order 
parameter - electric field i?3 and shear stress erg. Then the relevant order parameter, after which the 
phase transition takes place is a certain linear combination of P3 and eg- To determine it, a linear 
transformation {Pz,£q) — > (^,C) is performed, which diagonalizes the quadratic form (^3,66) in the 



expansion (2.1). Among ^ and C, the relevant order parameter is that one, which prefactor in (2.1) 
vanishes at the Curie point. 

However, in practical calculations for the only experimentally accessible field conjugate to the order 
parameter - the electric field E^, one usually uses the expansion in powers of polarization P3 

G = Go + °'(^ ~^°) p2 + ^ c , _ 
I 4 

and the higher order terms can be also taken into account. For the sake of simplicity, in this section we 
omit indices near polarization P3 and field Ej,. 

In the vicinity of the transition temperature at i? = 0, the thermodynamic potential G{P) has three 
minima P^^ = (— &± v6^-4ac)/2c and P2 = 0, the depths of which become equal at the transition point. 
Hence, the temperature of the first order phase transition Too at i5 = is found from the conditions 

r)G 

— = 0, G(Pi,3) = G(0), 



that are equivalent to 



3b^ — 16ac, Tqo = Tq 



362 



16ca' 

The extremum points of the thermodynamic potential G{P) aX E ^ Q are given by the equation 



E = aP + bP^ + cP"\ 



(2.3) 



Among five solutions of this equations, two correspond to maxima, and three Pi < 0, P2 ~ 0, and P3 > 
correspond to minima of G{P). From the condition 



G{P2) = G(P3) at £; > 0, G{Pi) = G(Pi) at £; < 0, 



in a Hnear in E approximation we find a field dependence of the first order phase transition temperature 
(the phase equilibrium curve) 

12c 

Tc = Tco — E- 



la'h 
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2.1 Critical point. Experimental measurements 

The phase equilibrium curves (the first order phase transition lines) terminate in two critical points, which 
coordinates (±i?*, T*) follow from the thermodynamic conditions - equation of state and the equilibrium 
condition 



dPjj^ dP dP^ 

In these points a difference between the truly ferroelectric and pseudo - "paraelectric" with a non-zero 
polarization P3 and strain eg phases disappears, so that the transition from one to the other takes place 
smoothly without jumps of the thermodynamic quantities. 

Experimentally the coordinates of the critical points are usually found using the phenomenological 
theory, by measuring in some way the coefficients of the Landau expansion. To determine the critical 
field E* and temperature T*, the following formulas are used 



E* 



a*P* + b{P*y + c{P*y + d{P*) 



T* = Tn 



To 



■ih{P*Y + bc{P*Y + 7d{P 



3* \6 



where the higher order terms in the expansion are also taken into account, and P* depends on the number 
of terms in the expansion 



0,^7^ 0; 




The simplest is to find the coefficient a', which is nothing but the inverse Curie constant of a crystal. 
The coefficients at higher terms in the expansion are determined by analyzing the field and temperature 
dependences of the crystal polarization. Usually, the analysis of the so-called F-plots is performed. If one 



restricts the Landau expansion by the P^ term, then from the equation for field (2.3) it follows |To|, that 
the dependence of the quantity 



r(p2 



E-aP 
P3 



b + cP^ 



on P^ must be linear and independent of the magnitude of the applied field. The intersection of this line 
with the ordinate axis (P^) gives a value of 6, and its slope gives the coefficient c. Hence, by analyzing 
the experimentally measured dependence r(P^) one is able to determine two coefficients of the Landau 
expansion. 

When the dependence r(P^) is closer to a quadratic one rather than to to linear one, we should 
neglect the term cP^ in the Landau expansion and take into account the term dP'', instead. Then the 
dependence 

r(p4)^:^^^5 + rfp4 



p3 



must be linear. 



Another, somewhat more sophisticated method of determination of the Landau coefficients was pro- 
posed by Schmidt Instead of the F-plots, the so-called isopols - lines of a constant polarization in 
the {E, T) plane should be analyzed. From (|2.3|) it follows that 



E bP^ + cP^ 

a'P a' 

that is, the dependence T{E) at constant polarization is linear. Approximating the experimental para- 
electric isopols, one can find the coefficient a' from the expression 



, 1 fdE\ 
a — — 



P \dT J 



p 



It should be noted that at the first order phase transition the isopols do intersect, whereas at the second 
order phase transition they do not. Therefore, this analysis permits one to determine the order of the 
phase transition and estimate the coordinates of the critical point (the point where the order of the phase 
transition changes), if one exists. 

Unfortunately, the described method of determination of the critical point coordinates according to 
the Landau coefficients does not allow to find more or less trustworthy values of the critical quantities. 
In Table |] we presented the measured in different papers Landau coefficients and calculated coordinates 
of the critical points for KH2PO4 (the Table is taken from Q). All the present results agree as for the 
magnitude of the critical field in this crystal about 100-300 V/cm. We did not present here the data 
obtained within dilatometric X-rays studies O], giving the critical field in KH2PO4 of about 8500V/cm. 



Table 1: Landau coefficients and calculated coordinates of the critical point for KH2PO4 (the data are 
taken from 0). 



Ref. 


a 


b 


c 


d 


E* 


T* - To 


method 




(lO-^esu) 


(lO-i^esu) 


(lO-i'^esu) 


(10-2^esu) 


(V/cm) 


(K) 




12 




3.9 


-1.9 


6.3 





120 


0.07 




13 




4.2 ±0.1 


-1.9±0.1 


5.4 ±0.4 





160 


0.07 


r-plot 


14 




4.3 ±0.2 


-2.35 ±0.4 


5.91 ± 1.5 





232 ± 70 


0.10 ±0.03 


Isopols 


15 




4.0 ±0.2 


-1.48 ±0.2 


3.1 ±0.4 





186 ± 60 


0.08 ± 0.03 


Isopols 




3.91 ±0.04 


-1.26 ±0.05 


3.2 ±0.1 





123 ± 18 


0.057 ± 0.007 


Isopols 




3.8 ±0.1 


-3.0 ±0.8 


6.5 ±1.1 





370 


0.16 


r-plot 




3.8 ±0.1 


-0.5 ±0.3 





3.8 ±0.4 


87 


0.036 


r-plot 




3.9 


-0.54 ±0.05 





2.85 ±0.10 


124 


0.046 






3.9 


-1.85 ±0.25 


3.3 ±0.5 


0.87 ±0.5 


280 


0.11 





As we see, a considerable discrepancy between values of the Landau coefficients takes place, and values 
of the critical point coordinates are very sensitive to the choice of the Landau coefficients. Thus, a 1.5 
times increase of the b coefficient leads to a 3 times as large value of the critical field and twice as large 
the quantity T* — Tq. The main sources of differences in the values of critical coordinates determined via 
the Landau coefhcients are the following 

1. Inadequacy of the Landau expansion as such. It is known to be suitable for the second order phase 
transitions, only. For crystals, where the first order phase transition is observed, a large number of 
terms in the expansion should be taken into account, and even then it is valid only in the narrow 
vicinity of the transition point. Even for KH2PO4, where the jump of polarization at the transition 
point is not as large as in deuterated KD2PO4, the Landau expansion is not quite adequate. 

2. Ambiguity of the Landau expansion. At interpretation of the same experimental data by different 
expansions (with c = 0, d 7^ or with c ^ 0, d = 0) the critical field can change to more than 4 
times. The expansion with c = 0, d 7^ yields considerably smaller values of the critical field than 
the expansion with 0, d — does. 



3. A large numerical error in calculations of the expansion coefficients at a given form of the expansion. 
Thus, at determining the coefficient b from F-plots, the error arises due to the data dispersion for 
r(P^) at small values of polarization. 

4. Changes in the Landau coefficients from sample to sample. First, a quite large dispersion of the data 
(as for determination of the critical point coordinates) takes place for the inverse Curie constant 
(the coefficient a') Similarly, the values of b for different samples of the crystal do differ, even 
though they are obtained by the same method. Here, the dispersion of the data may exceed the 
error of measurement |l^ . Apparently, the Landau coefficients and, thereby, the coordinates 
of the critical point, essentially depend on a quality of the sample, its history, presence of inner 
defects, etc. 

5. Different experimental methodics of polarization measurements (pyroelectric, dilatometric, hystere- 
sis loops) give values of the coefficient b, which may significantly differ. 

Comparison with experimental data of other quantities, calculated with the Landau expansion (for in- 
stance, temperature and field curves of polarization) does not allow to ascertain these coefficients, because 
these other quantities are not that sensitive to the values of the Landau coefficients, as the coordinates 
of the critical point. 

Hence, due to a large error in determination of the Landau coefficients and ambiguity of the Landau 
expansion, the method of calculation of the critical fields and temperatures via the phenomenologic theory 
allows one only to estimate the coordinates of the critical point (order of magnitude of the critical field). 



3 Microscopic model 

We consider a ferroelectric crystal of the KH2PO4 type to which external electric field E3, inducing 
polarization P3 and strain eg in the high-temperature phase, is applied. 

The entire Hamiltonian of the model |3j consists of a "seed" part, independent of a hydrogen subsystem 
configuration and attributed to a host lattice of heavy ions, and of pseudospin short-range and long-range 
hydrogen Hamiltonians, tunneling being neglected 



H ^ Nv 



Eo \ 

66 ^2 „0 p _ A33 p2 \ , TT 
~ ^36^3^6 ) + J^li 



hort • 



(3.4) 



The "seed" energy expressed in terms of the electric field E3 and strain eg includes the elastic, piezoelec- 
tric, and electric counterparts, c^q, Cgg, and X33 are the so-called "seed" elastic constant, coefficient of 
the piezoelectric stress, and dielectric susceptibility, respectively; v = v/kB, v is the primitive cell volume; 
fee is the Boltzmann constant; N is the number of primitive cells. 

-ffiong is the mean field Hamiltonian of the long-range dipole-dipole and lattice mediated p8| interac- 
tions between deuterons plus a linear in strain eg molecular field 0, 0| induced by piezoelectric coupling 



2 ^ ' 2 2 

q'f'qf 

2Nuif ~ {2vr^ - 27/>geg) 
9/ 
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E 
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(3.5) 



where 

Av = Jii + 2 J12 + Ji3 

is the eigenvalue of the long-range interaction matrix Fourier transform Jffi — ^ ^ Jff'{qq'); 



^q3) 



{'^qi) 



is the mean value of the Ising pseudospin a^f = ±1 which two eigenvalues are assigned to two equilibrium 
positions of a hydrogen on the /-th bond in the q-th unit cell. 



The Hamiltonian of the short-range configurational interactions between hydrogens is usually chosen 
such as to reproduce the energy levels of the Slater- type model for KDP (see, for instance Ref. ||l^) - the 
energy levels of up-down (twice degenerate at eg = and i?3 = 0), lateral (four- fold degenerate), 
single-ionized Si (eight-fold degenerate), and double- ionized Sq (twice degenerate) hydrogen configurations 
{Es < Ea < £i < eo)- 

Since the system is no longer symmetric with respect to the reflection (7^, in the ab plane and mirror 
rotation S4 around the c axis (both operations change the sings of polarization and strain), in presence of 
strain eg and in the electric field E3 a splitting of the energies of up and down {i = 1 and i = 2), lateral 
{i ~ 5, 6) and (i = 7, 8), and single ionized (i = 9, 10, 11, 12 and i — 13, 14, 15, 16) configurations takes 
place (confuguration numbers i are given in Appendix). Since strain eg and polarization P3 transform 
after the same irreducible representation {B2 in the paraelectric phase and Ai in the ferroelectric phase) , 
the field E3 does not split those levels which remain degenerate in presence of the strain eg. 

That part of the splitting, which origin is the distortion of the PO4 groups and changes of the angle 
between perpendicular in the paraelectric phase hydrogen bonds is described by linear functions of the 
strain eg. The contribution of the electric field E3 to the splitting is taken to be equal to the difference 
between the energy levels of the groups dipoles in the longitudinal field E3. Since the groups K-D„P04 
with hydrogens in up (i = 1) and down {i = 2) configurations and in single- ionized configurations 
(configurations i — 9, 10, 11, 12 and i — 13, 14, 15, 16) have oppositely directed projections of the dipole 
moments fi^ on the axis c, these groups have also different energies /XjE in external electric field . Usually, 
the absolute values of the c-projections of the dipole moments of up/down configurations /is are assumed 
to be twice as large as the corresponding projections of single-ionized configurations. The c-projections 
of dipole moments of lateral and double-ionized groups are zero. 

Finally, the Hamiltonian of the short-range interactions rewritten in terms of pseudospins according to 
standard rules |j and taken into account within the most appropriate for these crystals four-particle 
cluster approximation, is 
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(3.6) 



Detailed derivation of the Hamiltonian and the tabic of the split energy levels in presence of external 
electric field (longitudinal and transverse) are given in Appendix. 
In Eq. (3.6) the following notations are used 



y = [/ = - e, $ = 4e -I- 2wi - 8w; 



where 



£■ En e.s 



£1 



= eo - e^ 



are the so-called Slater energies. Splitting of the energy levels due to lowering the system symmetry is 
the only changes in the short-range interactions that we take into account within this model. Therefore, 
e, w, wi do not depend on strain eg of field E^. 

The condition of equality of the mean values of pseudospins rj = {<Jqf) calculated with the four-particle 
( ^ ) and single-particle deuteron Hamiltonians (see ||]) permits to exclude the self-consistency parameter 
A. The order parameter 77 and strain eg can be found by minimization of the thermodynamic potential 
(Gibbs' function) 5i£;(T, ag, S3, 77) 



5iij(T, ag, iJa, r,) = ^^^^4 - ve^e^E^, - -x^El 



2Tln2 + 2v'if - 2Tln(l - Tf )D - wgeg, 



where 



D = cosh(22; -I- /3(5sgeg) -I- 46cosh(2; — /3(5igeg) -f aag H h d, 

ag 



^ = o 1^ "i + /^^^ ~ /^^eee + — ^ — , 

2 1 — 77 2 

a = exp(— /?£), 6 = exp(— /3ti;), d = exp(— ag = exp{—(36ae£e)- 

The conjugate to the strain eq stress is introduced for the sake of simpUcity; in numerical calculations 
CT6 =0. 

Using the system of equations 



'-66 
„0 



£0 i 



4ipG m 2Mg 



-nj - eggee + X33-C'3 + 2— — , 
which follows from the thermodynamic equilibrium conditions 



3 / T,o 



V \ de& J T^E^^a, v\dE 
we find the isothermal dielectric susceptibility of a free crystal (at erg = const) 



Xl3 



^Ts , (635) 
>^33 + TE ' 

'-ee 



V5^3/T.e 

expressed via the isothermal dielectric susceptibility of a clamped crystal (at eq = const) 



X33 



the isothermal coefficient of piezoelectric stress 



^^36 



de, 



-6 / T,E3 

and the isothermal clastic constant at constant field 



,2^3 

^^36 ^ 



V D — 2mp 



(3.7) 



(3.8) 



(3.9) 



TE 
^66 



„i;o , S/jy-e + r , 2/3 2 

Cgg i ~ ~ r _ „n 



6 / T,_E3 



2^ 
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The other isothermal dielectric and piezoelectric characteristics can be recalculated via the found above 

by using the known formulas 

the constant of piezoelectric stress /136 



/i36 



dee, 



TP3 



\dPsJ 



Tee 



^36 
X33 



the elastic constant Cqq at constant polarization 



'-66 



( ^^f> \ 



C66 +636/136- 
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compliance at constant field 



TE 
'66 



the coefficient of the piezoelectric strain 
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daej 



TE3 



^TE' 
^66 
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(3.10) 



(3.11) 



(3.12) 



(3.13) 



the constant of piezoelectric strain 



536 = 



TPa 



'■36 

TP ■ 

66 



(3.14) 



When dynamic experimental methodics are used (measuring frequency exceeds the thermal relax- 
ation frequency), not the isothermal but adiabatic characteristics of the crystals are measured. In the 
paraelectric phase in absence of electric field the adiabatic and isothermal quantities coincide, but they 
may essentially differ in the ferroelectric phase or in the electric field. The adiabatic quantities can be 
obtained from the isothermal analogs with the help of the following formulas pO| 
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where 



S = R{~2 In 2 + 2 ln[l - r/^] + 2 In L> + 4r^^r; 



2M6 



) 



is the entropy of the hydrogen subsystem, p^ is the pyroelectric coefficient, cp^^a is the specific heat at 
constant polarization and pressure. From (3.15) it follows that the most prominent difference between the 
adiabatic and isothermal quantities takes place in the range of rapid temperature changes of polarization 
(where {dP^/ dT)\^ is large) - in the vicinity of the transition point or of the permittivity maximum. 

Using (3.15) and the relations (3.7)-(3.14) we can show that the adiabatic and isothermal "true" 
constants of the crystals - the elastic constant at constant polarization Cqq and piezomodules h^Q and g^Q 
- coincide. 

The pyroelectric coefficient and the contribution of the hydrogen subsystem into the specific heat at 
constant pressure and polarization are given by the following expressions 



Pe =P6 + 636016, 



and 
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(3.16) 
(3.17) 



where Acpa^ is the specific heat at constant polarization and strain 
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R is the gas constant; is the strain heat at constant P3 
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ae is the factor produced by temperature variation of the strain eg 
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(true thermal expansion, which invokes the diagonal strains, is not considered here. We used the following 
notations 



^6 = ^ 



(e + 5aQeQ)'^aae + (e - ^06^6)^ V ^bwdiQEQ smh.{z - I35iq£q)- 
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+4&f + (5jg£g] cosh(z — /J^igEe) + wi^d — {5s&£q)^ cosh(2z + P5a&£Q) 

96 = ^[-^s6£6 cosh(2z + (35as£&) + SMigEe cosh(z - pSiaee) + 26u;sinh(z - PSiqEq)] - ijMe, 
Ag = ^ [S'^qEq sinh(2z + /S^aeEe) + SfQEeib cosh{zQ - pSieEe}^ 

+Sa6 [aaeie + daeSe) - — (e - (^ae^e) ) + SiQweeAbsmh{z - (iSieee) 
\ ae J 

Mq — (3{Ahw cosh z + dwi + aeag + aea^^ + AfgEg). 
R is the gas constant. 

Having the pyroelectric coefficient and the specific heat at constant polarization, we can find the total 
molar specific heat at constant pressure 

A4-r(^||) - Acg^^ + ge'K- (3.18) 
, ( dS\ V 2RT ^ r ~ ^ 

is the polarization heat at given eg, 

4 Discussion 

4.1 Fitting procedure 

The presented in previous Section theoretical results will be used to description of the physical charac- 
teristics of a highly deuterated crystal KD2PO4. We shall also verify whether the presented theory is 
capable of describing behavior of the physical characteristics of an undeuterated KH2PO4 without taking 
into account tunneling effects. 

For a highly deuterated crystal K(Hi_2:Da;)2P04 with the transition temperature in zero field Tqq — 
211.73 K (a nominal deuteration level x = 0.89, hereafter abbreviated as KD2PO4) we use the values 
of the theory parameters found in |^ , providing a fair quantitative description of temperature of several 
associated with strain eg dielectric, piezoelectric, and elastic characteristics of the crystal at atmospheric 
pressure, namely, the dielectric permittivities of a free and clamped crystals, elastic constants and 

of open-circuited and short-circuited crystals, and piezomodules esg, dag, gag and /lag. Details of the 
fitting procedure are given in 

For and undeuterated crystal KH2PO4 we also chose the theory parameters such as to fit the tem- 
perature behavior of the elastic constants and c^^g, piezomodules eag, c?a6: 536, and /lag, as well as the 
longitudinal static dielectric permittivity of a clamped crystal ela. It should be noted that only the data 
of [2^ for the temperature dependences of the coefficient of piezoelectric strain dag , of H^, |6| 
for the dielectric permittiv ity eaa , and of p7| for the elastic constant are direct experimental points 



for KH2PO4. Using Eqns. ( |3.7D - ( |3.14D and having the values of dag, e'^^, and cfg, we find "experimental" 
points for the piezoelectric constants eae, /lae, 536, elastic characteristics c^^g, s^, s^q, and clamped dielec- 
tric permittivity Xaa- Recalculated thus values of the piezoelectric constants (/a6 ^ 53 • 10~^ cm^/esu and 
/lae ~ 3.8-10'' dyn/esu above the transition point are greater than 44-10"^ cm^/esu and 3.0 • 10* dyn/esu, 
respectively, given by Mason |28 but closer to the deQuervain values (50 • 10~^ cm^/esu for /laa)- 

By this set of the parameters we also try to describe the field dependences of polarization Pa 
longitudinal static dielectric susceptibility xt^ ^'^d elastic constant cf^ of a KH2PO4, obtained in It 
should be noted that in Q only the elastic constant was measured directly by the ultrasonic technique, 
whereas the susceptibility was recalculated using the known thermodynamic relations and the data of 
[ p9| for the piezomodules. 

A scheme of choosing the theory parameters for KII2PO4 was analogous to that for KD2PO4. At 
chosen e and w, close to those used in theories where the piezoelectricity and tunneling effects are not 



taken into account [|o|, and setting the Curie temperature of a clamped crystal JJ^'™? — 118.65 K, we 
found the value of the long-range interaction parameter v. Setting the values of the deformation potentials 
^6, <5a6j ^16, we found the value Ssq by fitting to experimental value of the piezomodule 635 at a certain 
temperature. The parameters e, w, Ve, SaQ, <5i6, c^^ , Ma/''^ were chosen such as to obtain the correct 
values of the transition temperature Tco = 122.8 K Q of a free crystal, ratio of the polarization jump at 
the transition point to saturation polarization (jump of the order parameter), correct temperature curves 
of the elastic constants [p7| a nd c£, piezomodules c?36, 536, and h^Q, longitudinal static dielectric 
permittivities of a free efg p5|, |37[ |2^ , gj, Q and clamped £^3 crystals, as well as specific heat. Properly 
chosen parameters egg, Xgg'' yield the correct values of the piezomodules d^Q, egg, (73g, and /136 in the 
high-temperature limit. 

In Table || we present the adopted values of the theory parameters for KD2PO4 and KH2PO4. 



Table 2: The theory parameters. 





e w V tpe 5s6 Sae (5i6 
(K) 


Cgg ■ iU 636 X33 

V 

(dyn/cm^) (/iC/cm^) (esu/cm^) 


KD2P04 
KH2P04 


91.3 781 34.615 -500 -692.5 1350 50 
52.5 383 20.58 -380 -257 300 40 


6.7 6.9 0.42 • 10" 0.4 
7.3 6.4 0.3 • W 



4.2 Effective dipole moment 

In our calculations, as well as in earlier theories (see, for instance, 1^) two different values of the effective 
dipole moment /Ltg < must be used in the ferroelectric and paraelectric phases. The value of 2/^*3/1; 
in the ferroelectric phase is given by the experimental value of saturation polarization. The paraelectric 
values of 2/xj^/z;, listed in Table ^ are set by the relation between the experimental and theoretical values 
of the Curie constant for the longitudinal static dielectric permittivity. In the present theory, /itg must 
also provide an agreement with experiment for the piezomodules 

636 A^3 7 d^Q ^ ^3, ^36 ~ ^, 936 ^ ^ 

in the paraelectric phase. In KH2PO4 and KD2PO4 the ratio is about 1.2. Neither in the present 

model nor in simpler versions of the proton ordering models it is possible to choose the theory parameters 
such as to fit the temperature curves of polarization and permittivity with a single value of the parameter 
M3- 

The necessity to use in calculations two different values of jj, for paraelectric and ferroelectric phases 
is not a problem while we do not consider behavior of the system in external fields conjugate to the 
order parameter. These fields smear out the phase transition; therefore, when the difference between the 
phases becomes only quantitative (in fields above the critical one), it does not make sense to distinguish 
the paraelectric and ferroelectric phase and the corresponding /x"*" and /i^. 

Tokunaga |]3l| relates the deviation of the ratio /x"*"//^" from unity with the existence of an un- 
derdamped soft mode in a crystal. In the displacive type ferroelectrics ^ 1, whereas in the 

order-disorder type ferroelectrics fi~^ / fi'' ^ 1. Failure of the order-disorder model to describe both po- 
larization and dielectric permittivity without invoking two values of /x for different phases shows the 
limits of the model suitability. Since the KH2PO4 family crystals undergo phase transitions of a mixed 
ordering-displacive type (it is believed that hydrogen ordering triggers displacements of heavy ions), for 
a consistent description of dielectric properties of these crystals, we need to sophisticate the model and 
include phonon degrees of freedom and anharmonicity into consideration. Moreover, we need to go fur- 
ther than Kobayashi model (mixed proton-phonon model with only one optical mode taken into account) , 
which is known only to renormalize parameters of the proton ordering model. 



4.3 Zero field case 



Figures H and || illustrate how the proposed theory describes experimental data for the physical charac- 
teristics of KD2PO4 and KH2PO4 crystals in absence of external electric field. Polarization of a pure 
KH2PO4 P3 is calculated with 2fi^ /v = 5.0 /iC/cm. To obtain the total specific heat, we added a con- 
stant term 60 J/mol K, describing the background specific heat of a host lattice of heavy ions, to the 



calculated from (3.f8) contribution of the hydrogen subsystem to the specific heat. 



4.4 Non-zero field case 

In figure || we plot the Tc — phase diagrams of the KD2PO4 and KH2PO4 crystals. These diagrams 
are of the same topology as the Tq — ce diagram of KD2PO4 The coordinates of the critical point, 
calculated within the microscopic theory, essentially depend on the magnitude of the order parameter 
jump Ar] (the ratio of the polarization jump to the saturation polarization) at the transition point in 
zero field. Since the order of the phase transition in a pure KH2PO4 is close to the second one, while in 
a deuterated KD2PO4 a pronounced jump of polarization at the transition point takes place, the critical 
field in KH2PO4 is much lower than that in KD2PO4. 

For KD2PO4 we obtain a good agreement with experimental data for a field dependence of the 
transition temperature dTc/dE^ — 0.13K/ kVcm in our model and 0.125 ± 0.1 in [^. Theoretical 
value of the critical field — 7.0 kV/cm, where the phase transition in the system disappears, accords 
fairly well with the experimental estimate 7.1 ± 0.6 kV/cm The critical temperature T* = 212.6 K is 
close to that from the Tq — diagram. 

For an undeuterated crystal, the theory gives the following characteristics of the phase diagram 

E* = 105 V/cm, T* -To = 0.05 K, dTc/dE^ = 0.23K/ kVcm. 
The calculated critical field in KH2PO4 agrees with the the estimates 80 ^ 370 V/cm obtained via 
experimentally measured coefficients of Landau expansions. However, at varying the microscopic Arj 
within experimentally estimated limits (for which slight changes in the theory parameters are required), 
the calculated critical field is significantly increased, while more or less precisely measured macroscopic 
characteristics (e.g. permittivity) remain practically unchanged. We may conclude that the microscopic 
theory is able only to estimate the coordinates of the critical point. 

In addition, in fields close to critical, peaks of the longitudinal permittivity or other characteristics 
having peculiarities at the transition point are very sharp. Therefore, even small changes of temper- 
ature (-- 10-3 K) can change the permittivity value several times. The maximal temperature step in 
calculations, yielding more or less correct value of the permittivity maximum, is ^ 10~* K. 

As we have already mentioned, we cannot describe polarization of the crystals, without introducing 
the "ferroelectric" values of the dipole moment 113. Nevertheless, we can compare the field behavior 
of the order parameter rj with the experimental data for the ratio P^/Ps- Here Pg ~ 5.0 /iC/cm^ for 
KH2PO4 and 6.2 /iC/cm^ for KD2PO4 are the experimental values of saturation polarization, practically 
independent of external field [||. Such comparison is justified, since the contributions to polarization from 
the terms X33E3 &n.d e^gEg are not essential. Figure ^ shows a fair accordance of theoretical temperature 
and field dependences of Ps/Ps for KH2PO4 and KD2PO4 with the experimental data. Increase of 
polarization and vanishing of polarization jump are well reproduced by the theory. 

In figure ^ we depicted the temperature dependences of two characteristics - isothermal static di- 
electric permittivity of a free crystal e^^ and the piezomodule h^Q of a deuterated crystal at different 
values of external electric field. Temperature curves of other characteristics of the crystal exhibiting 
peculiarities at the transition point - compliance and piezomodules d^^ and 635 - are similar to those 
of the permittivity e^^. Temperature and field behavior of the piezomodule g^e to those of h^g. The 
corresponding curves for a undeuterated are qualitatively similar, however, due to closeness of the order 
of the phase transition to the second one, the temperature and field intervals where all these changes take 
place are much more narrow. 

Maximal values of e"^^ and similar quantities shift to higher temperatures with field. As the field 
approaches its critical value, the peak values of these quantities in the transition points increase and are 
maximal at the critical point. The higher fields smear out the phase transition and round off and lower 
down the peaks. 
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Figure 1: Temperature dependences of the strain eg-related physical characteristics. Lines and solid 
symbols correspond to a deuterated KD2PO4; open symbols representing pure KH2PO4 are shown for 
comparison. Experimental points are taken from □ - [|8|; o - j2^; o - ||2l|, |2^; ■ - [52| ; • - |33 ; A and 
▲ are recalculated from Eqns. ( |3.7| ) - ( ^.14 ) using experimental data of Refs. [^8[ and Refs. 
respectively. 
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Figure 2: Temperature dependences of the physical characteristics of an undeuterated crystal KH 
Experimental points: for polarization: □ - o - [Q; A - [^; V - iQ; for the permittivity: □ 
o - H; A - ||; V - H; o - |§; for specific heat: □ - ||; o - |3|; o - ||; for the others: □ 
o - |27|; V are recalculated from Eqns. (3.7) - (3.14) using experimental data of Refs. 
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Figure 3: Tc — -E3 phase diagrams of KD2PO4 and KH2PO4. Experimental points are talcen from 
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Figure 4: Ratio of polarization to saturation polarization P^/Ps as a function of temperature for KD2PO4 (left) 
and KH2PO4 (right) at different values of electric field E3 (kV/cm) (left): 1 - 0; 2 - 2.82; 3 - 5.64; 4 - 8.46; 5 
- 11.28: (right): 1 - 0; 2 - 5.813; 3 - 12.5 4 - 20.31. Experimental points are taken from for KH2PO4 and 
from lli for KH2PO4. 



The field effects on piezomodules h^Q and g^Q is similar to its effects on polarization. On increasing 
the field, jumps of these quantities at the transition point decrease and vanish at the critical field. At 
higher fields, these quantities exhibit smooth temperature dependences. 

It is interesting to note that field dependences of adiabatic quantities are stronger that those of the 
isothermal quantities. Thus, as fig. ^where we present temperature curves of the adiabatic and isothermal 
piezomodule dfg of a pure KH2PO4 shows, the peak values of the adiabatic piezomodule (as well as of 
Sgff , efg, £33) decrease with the external field (above the critical point) much stronger than the peak 
values of the isothermal d|"g, and the maxima of the adiabatic quantities shift with the field much more 
perceptibly that those of the isothermal quantities. 

As one can see (Fig. |^), the proposed theory yields a fair agreement with experiment for the field 
effects on the adiabatic longitudinal dielectric susceptibility of a free crystal and elastic constant cfff 
in a pure KH2PO4 0. Deviation of the theoretical curves from experimental points in the vicinity of 
the "transition" (maximum of the permittivity) is attributed to our neglecting tunneling effects. Such 
neglecting worsens an agreement with an experiment for the temperature curve of order parameter P3 / Pg 
and, thereby, of xfs and cf^. 
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Figure 5: Isothermal static dielectric permittivity a free crystal £33" and isothermal piezomodule /lae for a 
deuterated KD2PO4 as functions of temperature at different values of electric field E3 (kV/cm): 1 - 0; 2 - 3.82; 
3 - 5.64; 4- E* = 7.0; 5 - 8.46; 6 - 11.28 . 
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Figure 6: Isothermal and adiabatic piozomodules dse of a undeutcrated KH2PO4 as functions of temperature at 
different values of electric field E3 (kV/cm): 1 - 0; 2 - 1; 3 - 2; 4 - 3; 4 - 5. 

5 Concluding remarks 

In this paper we apply the previously developed microscopic theory to description of effects produced by 
the longitudinal electric field E3 on the phase transition and physical properties of the KH2PO4 family 
ferroelectrics. Calculations arc performed in the four-particle cluster approximation within the proton 
ordering model without taking into account the tunneling effects. The developed model takes into account 
the existence of the shear strain eq, which is spontaneous in the ordered phase and induced by the stress 
(76 and electric field (via the piezoeffect) in the disordered phase. We also take into account other 
effects induced by this strain - the splitting of short-range proton correlation energies and appearance of 
effective molecular fields, created by piezoelectric coupling. 

The presented model provides a fair agreement with experimental data for temperature behavior 
(at zero field) of the dielectric, piezoelectric, elastic, and thermal characteristics - second derivatives of 
thermodynamic potentials (dielectric permittivity, piezomodules, elastic constants, specific heat) of the 
crystals. Such agreement takes place for a deuterated crystal KD2PO4 and for a pure KH2PO4 without 
taking into account the tunneling. 

The model satisfactorily reproduced the parameters of experimental Tc — E3 phase diagrams of 
KD2PO4 and KH2P04. The transition temperature increases with field; this is accompanied by increase 
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Figure 7: Adiabatic elastic constant cffp and the inverse longitudinal static dielectric susceptibility of a free 
crystal (xfa')"^ for a- pure KH[2P04 as functions of temperature at different values of electric field E3 (kV/cm). 
Experimental points are taken from H]. 



in polarization in the high-temperature phase and decrease of its jump at the transition point. The 
phase equihbrium curves terminate in the critical points, where the polarization (order parameter) jump 
vanishes. At higher fields the temperature curves of polarization and strain eg are smooth. Such behavior 
is typical for the systems undergoing the first order phase transitions in external fields conjugate to the 
order parameter. 

The magnitude of the critical field essentially depends on the degree of the "first-orderness" of the 
phase transition). Thus, in KH2PO4 where the phase transition is close to the tricritical point, the 
critical field is several times smaller than in KD2PO4, where a pronounced first order phase transition 
takes place. Experimental method of the critical coordinates determination based in the phenomenologic 
analysis allows only to estimate the critical field. The calculated within the microscopic theory value 
of the critical field is determined by the magnitude of the order parameter at the transition point. 
Since on changing this quantity (at the proper and slight change of the relations between the theory 
parameters) within the limits of experimental estimates the calculated critical field significantly increases, 
the microscopic approach also permits only to estimate the coordinates of the critical point, though this 
estimate is better than the phenomenologic one. 

As the field approaches the critical point, the peak values of thermodynamic quantities that have 
peculiarities at the transition point (the longitudinal dielectric permittivity £33, piezomodules d^Q and 635, 
and the elastic compliance Sqq) increase, while the jumps of "true" constants of the crystals (piezomodules 
536 and ft,36, and elastic constant Cq^) decrease. At the critical point the peak values of the formers 
are maximal, whereas the jumps of the latters vanish. In fields above the critical one, temperature 
dependences of all thermodynamic quantities are smooth, and the maximal values of the characteristics 
that have peculiarities at the transition point decrease. 

In presence of external field, there arise a difference between isothermal and adiabatic characteristics 
(dynamic methods of measurements yield the adiabatic quantities) , which have peculiarities at the tran- 
sition point. The field dependences of the adiabatic quantities are stronger than those of the isothermal 
quantities. Thus, in fields above the critical one, the peak values of the adiabatic quantities decrease 
with the field much faster that the peaks of the isothermal quantities, and the maxima of the adiabatic 
quantities shift to higher temperatures that the maxima of the isothermal quantities at the same field. 
Within the proposed theory we obtain a fair description of the experimental data for the field dependences 
of the adiabatic dielectric susceptibility and elastic constant of KII2PO4. 

The developed model, as well all earlier used versions of proton ordering model, requires two different 
values of effective dipole moment to be introduced for description of polarization and dielectric permit- 
tivity of the crystals (in paraelectric and ferroelectric phases). For a consistent description of dielectric 
properties of these crystals, phonon degrees of freedom and anharmonicity must be taken into account. 
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